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Abstract Large-scale structure formation, accretion and merging processes, AGN ac- 
tivity produce cosmological gas shocks. The shocks convert a fraction of the energy 
of gravitationally accelerated flows to internal energy of the gas. Being the main gas- 
heating agent, cosmological shocks could amplify magnetic fields and accelerate ener- 
getic particles via the multi-fluid plasma relaxation processes. We first discuss the basic 
properties of standard single-fluid shocks. Cosmological plasma shocks are expected to 
be collisionless. We then review the plasma processes responsible for the microscopic 
structure of collisionless shocks. A tiny fraction of the particles crossing the shock is in- 
jected into the non-thermal energetic component that could get a substantial part of the 
ram pressure power dissipated at the shock. The energetic particles penetrate deep into 
the shock upstream producing an extended shock precursor. Scaling relations for post- 
shock ion temperature and entropy as functions of shock velocity in strong collisionless 
multi-fluid shocks are discussed. We show that the multi-fluid nature of collisionless 
shocks results in excessive gas compression, energetic particle acceleration, precursor 
gas heating, magnetic held ampliflcation and non-thermal emission. Multi-fluid shocks 
provide a reduced gas entropy production and could also modify the observable ther- 
modynamic scaling relations for clusters of galaxies. 
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1 Introduction 



The observed large scale structure of the Universe is thought to be due to the grav- 
itational growth of density fluctuations in the post-inflation era. In this model, the 
evolving cosmic web is governed by non-linear gravitational growth of the initially 
weak density fluctuations in the dark energy dominated cosmology. The web is traced 
by a tiny fraction of luminous baryonic matter. Cosmological shock waves are an essen- 
tial and often the only way to power the luminous matter by converting a fraction of 
gravitational power to thermal and non-thermal emissions of baryonic/leptonic matter. 

At high redshifts {z > flOO) the pre-galactic medium was hot, relatively dense, 
ionised, with a substantial pressure of radiation. The cosmic microwave background 
(CMB) observations constrain the amplitudes of density inhomogeneities to be very 
small at the last scattering redshift z ~ 1000. Strong non-linear shocks are therefore 
unlikely at that stage. The universe expands, the matter cools, and eventually recom- 
bines, being mostly in neutral phase during the "dark ages" of the universe. At some 
redshift, 6 < z < 14, hydrogen in the universe is reionised, likely due to UV radia- 
tion from the first luminous object s, leaving the intergalactic medium (IGM) highly 
reionised (see e.g. iFan et al] (|2006| ) for a recent review). The reionisation indicates 
the formation of the first luminous objects at the end of the "dark ages", either star- 
forming galaxies or Active Galactic Nuclei (AGN). The compact luminous objects with 
an enormous energy release would have launched strong (in some cases, relativistic) 
shock waves in the local vicinity of the energetic sources. At the same evolution stage, 
formation of strong density inhomogeneities in the cosmic structure occurs. Since then 
the non-linear dynamical flows in the vicinity of density inhomogeneities would have 
created large scale cosmic structure shocks of modest strength, thus heating the bary- 
onic matter and simultaneously producing highly non-equilibrium energetic particle 
distributions, magnetic fields and electromagnetic emission. 

Most of the diffuse X-ray emitting matter was likely heated by cosmological shocks 
of different scales. Accretion and merging processes produce large-scale gas shocks. 
Simulations of structure formation in the Universe predict that in the present epoch 
about 40 % of the normal baryonic matt er is in the Warm-Hot Intergalactic M edium 
(WHIM) at overdensities 5 ~ 5 - 10 (e.g. lCen fc Qstrikei|[l999l : lOave et allboOll ). The 

WHIM is likely shock-heated to temperatures of 10^ 10^ K during the continuous 

non-linear structure evolution and star-formation processes. 

The statistics of cosmological shocks in the large-scale structure of the Universe 
were simulated in the cont ext of the TlCDM-cosmology using PM / Eulerian adia batic 
hydrodynamic codes (e.g. iMiniati et alllioool : IRvu et all 12003): iKang et al.lboOTi') an d 



more recently with a smoothed particle hydrodynamic code by Pfrommer et al. ( 2006l l 



They identified two main populations of cosmological shocks: (i) high Mach number 
"external" shocks due to accretion of cold gas on gravitationally attracting nodes, and 
(m) moderate Mach number (2 < A4s < 4) "internal" shocks. The shocks are due 
to supersonic flows induced by relaxing dark matter su bstructures in re latively hot, 
already shocked, gas. The internal shocks were found bv lRvu et all ^2QQi ) to be most 
important in ener gy dissipation providing intercluster medium (ICM) heating, and they 
were sugg ested bv lBvkov et"all (|200d ) to be the likely sources of non-thermal emission 
in clusters of galaxies. 

Hydrodynamical codes deal with N-body CDM and single-fluid gas dynamics. How- 
ever, if a strong accretion shock is multi-fluid, providing reduced post-shock ion tem- 
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perature and entropy, then the internal shocks could have systematically higher Mach 
numbers. 

Space plasma shocks are expected to be collisionless. Cosmological shocks, being 
the main gas-heating agent, generate turbulent magnetic fields and accelerate energetic 
particles via collisionless multi-fluid plasma relaxation processes thus producing non- 
thermal components. The presence of these non-thermal c omponents may affect the 
global dynamics of cluster s of ga laxies [Ostriker et al.l (|2005| ) and the (Tv-T, M-T, Lx- 
T scaling relations iBvkovl (l2005h . Detailed discussion of the cosmological simulations 
of the scaling relat ions with account of only thermal components can be found in 
iBorgani et al.l[2008l - Chapter 13, this volume. 

In Sect. [2] we discuss the basic features of the standard coUisional shocks. The main 
part of the review is devoted to physical properties of cosmological shocks with an 
accent on collisionless shocks and associated non-thermal components. In Sect. [5] we 
discuss the most important features of multi-fluid collisionless shocks in the cosmo- 
logical context including the effects of reduced entropy production, energetic particle 
acceleration and magnetic field amplification in the shocks. 



2 Single-fluid MHD-shocks 



Shock waves are usually considered as a sharp transition between a macroscopic su- 
personic (and super- Alfvenic) upstream fiow (state 1) and slowed down to a subsonic 
velocity downstream fiow (state 2), providing a mass flow jn through the shock sur- 
face. It is assumed that a gas particle (or an elementary macroscopic fluid cell) is at 
any instant of time in the local thermodynamic equilibrium state corresponding to the 
instantaneous values of the macroscopic parameters. The Maxwellian distribution of 
all species is ensured after a few molecular (or Coulomb) collisions have occurred. The 
macroscopic parameters characterising the state of the gas, such as density, specific 
internal energy, or temperature, change slowly in comparison with the rates of the 
relaxation processes leading to thermodynamic equilibrium. We consider here a single- 
fiuid plasma model assuming complete electron-ion relaxation. Under these conditions, 
in a frame moving with the shock front, with the matter fiux across the shock surface 
jn 7^0, the conservation laws for mass (in non-relativistic fiows) , momentum and energy 
can be written as follows: 



jn [^] = Bn [Ut] , (1) 

[-+P+^]-0, (3) 

+4 +P- -^^tUt]-0. (4) 

2p^ 2 47rp 47rjn 

Here U — {un, ut) is the bulk velocity, w — e + P/p is the gas enthalpy, e ,P,p 
are the internal gas energy, pressure and density respectively. The subscripts n and t 
are used for the normal and transverse components respectively. We used the standard 
notations [^1] — A2 ~ Ai for the jump of a function A between the downstream and 
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MHD- shock geometry 
(the coplanarity theorem) 



Fig. 1 A sketch illustrating the coplanarity theorem for a plane ideal MHD-shock. The up- 
stream and downstream bulk velocities Ui and U2, magnetic fields Bi and B2 and the shock 
normal N all lie in the same plane. The shock is at rest in the reference frame where also uti 
= 0. The shock is of infinitesimal width in the sketch. Simulated structure of the transition 
region of a coUisionless shock is shown in Fig. [2] and Fig. |3] where its finite width is apparent. 



up stream regions. In the MHD c ase the relations equivalent to Eq. [T]-|3] were obtained 
bv Ide Hoffmann fc TeUerl (|l95(i ). The equations are valid in the case of a magnetic 



field frozen-into moving plasma, where E = — ^ x B. A specific feature of MHD shock 
waves is the so-called coplanarity theorem (e.g. iLandau fc Lifshit z 198^) saying that 
the upstream and downstream magnetic fields Bi and B2 and the shock normal all lie 
in the same plane as it is illustrated in Fig. [T] It is important to note that if Bn 7^ 
there is an especial reference frame where local velocity U and magnetic field B are 
parallel both in the upstream and downstream, providing E = 0. 

From Eq. [l]-|4]one may obtain a generalised Rankine-Hugoniot (RH) adiabat 

£2 - ei + - + Pi) + ^{Bt2 " Bt^f)} = 0. (5) 

2 p2 pi 87r 

The RH adiabat connects the macroscopic parameters downstream of the flow once 
the upstream state is known. In a parallel shock (Bf=0) 

P2^ (7g + l)Ml 
PI (7g-l)M| + 2' 

72 ^ [27g.M'-(7g-l)][(7g-l).Ms+2] 
Ti (7g + l)2Mi ' 
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where 7g is the gas adiabatic exponent. We restrict ourselves here to a fast mode shock 
where Cgi < mi, and Da2 < U2 < Cs2, for Wai < Cgi- The phase velocity Via is the 
Alfven velocity in the downstream, Cgi, Cs2 are the sound speeds in the upstream and 
downstream respectively. We define here the shock Mach numbers as A^s = '^sh/csl 
and Ma. = iigh/^al- 

In the case of a perpendicular shock [Bn = 0) the compression ratio is 

r^Pl = + 

PI V + + 4(7g + 1) (2 - 7g) Ma ' ^ 

In a single fluid strong shock with A^s 3> 1 and Ma 3> 1 one gets 

T2~2- P^^-^ = 1-38 X 10^ vis (K), (9) 

(7g + 

for any magnetic field inclination (e.g. lOraine fc McKee|[l993l ). The mass per particle 
^ was assumed to be [1.4/2.3]mjj and Hgg is the shock velocity in 10* cms^^. 

The RH adiabat does not depend on the exact nature of the dissipation mechanisms 
that provide the transition between the states 1 and 2. It assumes a single-fluid motion 
in regular electromagnetic fields. However, the dissipative effects control the thickness 
of the shock transition layer. In the case of a weak shock of Mach number Ms — 
the thickness is large enough, allowing a ma croscopic hyd rodynamical description of 
the fluid inside the shock transition layer fe.g. lLandau fc L ifsliitz 1984|). The gas shock 
width A in collisional hydrodynamics without magnetic fields is given by 



" {P2-Pl){d'^V/dP-^)s' 

Following iLandau fc Lifshit3 (|l959l ) the gas shock width in Eq. [TD] can be expressed 
through the viscosities r] and C,, and thermal conductivity k, since 

14 11 

:[(^'? + C) + ^(---)] 



2pvi 3 Cv Cp 

Here Cv and Cp are specific heats at constant volume and at constant pressure respec- 
tively. Extrapolating Eq.[lO]to a shock of finite strength where P2 — Pi ^ P2, one may 
show that the gas shock width A is of the order of the mean free path A. 

It is instructive to note that the entropy is non-monotonic inside the finite width 
of a weak gas shock (A^s — 1) ^ 1 and the total RH jump of the entropy As across 
the shock is of the third order in {Ms — 1): 

1 /aV 



(1^) {P2~Plf <x{Ms^lf, (11) 



i2ri KdP'^J 



while the density, temperature and pressure jumps are oc {Ms — 1) ( Landau fc Lifshit3 
\l95i ). 

In plasma shocks the shock structure is more complex because of a relatively slow 
electron-ion temperature relaxation. Such a shock consists of an ion viscous jump and 
an electron-ion thermal relaxation zone. In the case of plasma shocks the structure 
of the ion viscous jump is similar to the single fluid shock width structure discussed 
above and can be studied accounting for the entropy of an isothermal electron fluid. 
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The shock ion viscous jump has a width of the order of the ion mean free path. The 
scattering length (the mean free path to 7r/2 deflection) Ap of a proton of velocity vj 
(measured in 100 kms~^) due to binary Coulomb collisions with plasma pr otons of 
density n (measured in cm~'^) can be estimated as Ap « 7 x 10^** v^n~^ cm ( Spitzej 
1 19621 ). After the reionisation [z < 6) the Coulomb mean free path in the WHIM of 
overdensity 5 is Ap 3.5 x 10^^ S^^ (1 + zy^^O^.h^ /0.02)^^ cm. Here and below 
i7l3 is the baryon density parameter. The mean free path due to Coulomb collisions is 
typically some orders of magnitude smaller than that for the charge-exchange collisions 
in the WHIM after reionisation. The ion-electron thermal relaxation occurs on scales 
about Ac X -^Z nip/Zmc. Since Ac ~ Ap, the width of the relaxation zone is substantially 
larger than the scale size of the ion viscous jump. The application of the single fluid 
shock model Eq. [TJ-jl] to electron-ion plasmas assumes full ion-electron temperature 
relaxation over th e shock width. For a discussion of the relaxation processes see e.g. 
iBvkov et aLlbOoil - Chapter 8, this volume, and references therein. 

In a rarefied hot cosmic plasma the Coulomb collisions are not sufficient to pro- 
vide the viscous dissipation of the incoming flow, and collective effects due to the 
plasma flow instabilities play a major role, providing the coUisionless shocks, as it 
is di rectly observed in th e hcliosphere. The observed structure of supernova remnants 
fe.g. lWeisskopf &; Hughe si i2006 l is consistent with that expected if their forward shocks 
are coUisionless. Moreover, the non-thermal synchrotron emission seen in radio and X- 
rays is rather a strong argument for high energy particle acceleration by the shock that 
definitely favours its coUisionless nature. That allows us to suggest that cosmological 
shocks in a rarefied highly ionised plasma (after the reionisation epoch) are likely to be 
coUisionless. There are yet ve ry few observational studies of cosmological shocks (e.g. 
iMarkevitch fc VikhUninlboOTi ). We review some basic principles of coUisionless shock 
physics in the next section. 



3 CoUisionless shocks 

Since the discovery of the solar wind in the early 1960's it has been realised that the 
rapid rise time of magnetic storms observed in the Earth suggested v ery thin coUi- 
sionl ess shocks created by solar flares (see for a discussion Sa gdeev..l96S : iKennel et al.l 
Il985l '). The thickness of a viscous jump in a strong collisional shock is of the order of a 
mean free path (see e.g. IZel'dovich fc Raizerlll967l ). The Coulomb collision mean free 
path in the tenuous solar wind plasma is comparable to the Sun-Earth distance, and 
thus the magnetic storm rising time due to standard collisional shocks would exceed 

the observed time by orders of magnitude. 

There are very specific features of coUisionless plasma shocks I SagdeevI Il966l ). 



Shocks in dense enough plasma with frequent Coulomb collisions evolve very fast to 
Maxwellian particle distributions with very few particles at high energies. On the con- 
trary, in coUisionless plasma shocks, a small minority of particles could gain a dis- 
proportionate share of the energy and become non-Maxwellian. CoUisionless shocks 
enable acceleration of a small fraction of the particles to very high energies. Moreover, 
the accelerated particles could carry away a substantial amount of the kinetic energy 
of the plasma flow dissipated at the shock. The energetic particles can penetrate far 
into the shock upstream gas, to create an extended shock precursor. The cold gas in 
the shock upstream is decelerated and pre-heated by the energetic particle and MHD- 
wave pressure on a scale larger than a mean free path of an energetic particle. This 
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occurs not only at the bow shock of the Earth at mode rately low en eraies, but also in 
astrophysical shocks at highly relativistic energies fe.g. lRussellll2005l l. 

A direct study of coUisionless shock waves in a laboratory is an extremely difficult 
task. Most of the experimental data on coUisionless shock physics are coming from 
space experiments. There are direct observational data on the shock wave structure in 
the interpla netary medium with clear evidence fo r ion and electron acceleration by the 
shocks (e.g. iTsurutani fc Linlll985l : iRusselllboOsI l . 

Computer simulations of the full structure of coUisionless shock waves describe 
the kinetics of multi-species particle flows and magneto-hydrodynamic (MHD) waves 
in the strongly-coupled system. The problem is multi-scale. It requires a simultaneous 
treatment of both "microscopic" structure of the subshock at the thermal ion gyroradii 
scale where the injection process is thought to occur, and an extended "macroscopic" 
shock precursor due to energetic particles. The precursor scale is typically more than 
10® times the microscopic scale of the subshock transition region. 

Energetic particles could be an essential component in the WHIM and clusters 
of galaxies. Nonthermal particle acceleration at shocks is expected to be an efficient 
process at different evolutional stages of clusters. Being the governing process of the 
supernova remnant coUisionless shock formation, nonlinear wave-particle interactions 
are responsible for both shock heating and compression of the thermal gas, as well as 
for creation of an energetic particle population. 



3.1 Micro processes in coUisionless shocks 

In the strong enough coUisionless shocks (typically with a Mach number above a few) 
resistivity can not provide energy dissipation fast enough to create a standard shock 
transition (e.g. iKennel et al]|l985l ) on a microscopic scale. Ion instabilities are impor- 
tant in such shocks that are called supercritical. 

At the microscopic scale the front of a supercritical shock wave is a transition region 
occupied by magnetic field fluctuations of an amplitude SB/B ~ 1 and characteristic 
frequencies of about the ion gyro-frequency. Generation of the fluctuations is due to 
instabilities in the interpenetrating multi-flow ion movements. The width of the tran- 
sition region of a quasi-parallel shock wave reaches a few hundred ion inertial lengths 
defined as 1^ — c/ajpj « 2.3 x 10^ n~^'^ cm. Here is the ion plasma frequency and n 
is the ionised ambient gas number density measured in cm~^. The ion inertial length in 
the WHIM can be estimated as k « 5.1 x 10^" 5~^^^ (1 + z)'^/"^ {Oi.h'^ /0.02)~^^^ cm, 
providing the width of the coUisionless shock transition region is smaller by many orders 
of magnitude than the Coulomb mean free path (that is in the kiloparsec range) . 

The transition region of a quasi-perpendicular shock is somewhat narrower. The 
wave generation effects at the microscopic scale have been studied in some detail with 
hybrid code simulations (e.g. lQuestlll988h . The large-amplitude magnetic field fluc- 
tuations in the s hock transition region were directly measured in the interplanetary 
medium (see e.g. iKan et al.|[T99ll ). 

There are a few ways to simulate numerically the kinetics of the coUisionless plasma 
phenomena. The most comprehensive study of the coUisionless shock structure can be 
performed with the particle-in-cell (PIC) method where all the plasma components 
are considered as discrete particles in self-consistent fields. The PIC method allows 
one to resolve electron scale lengths and frequencies, but on the other hand it requires 
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considerable computer resource. A serious constraint on PIC and other plasma parti- 
cle simulations of collisionless shock s is that they must be done fully in three spatial 
dimensions (3D). Ijones et al.l ( 1998h have proved that PIC simulations with one or 
more ignorable dimensions artificially confine particles to field lines and particularly 
eliminate cross-field diffusion. The effect is especially important for simulations of a 
creation of a superthermal particle population. All three box dimensions must be in- 
volved in these simulations. Exact modelling of electron kinetics in colli sionless shocks 
requne PIC simulations (e.g. .Hoshino fc Shimada ,20o3 : ISchmitz et al.|[2002i '). On the 
other hand, the bulk of the energy of a collisionless shock is carried by the ions and ve- 
locity relaxation processes are typically longer than the ion gyro-periods. Thus, though 
the basic shock physics evolve on ion spatial and temporal scales, the electron kinetic 
description requires fine resolution at electron scales. 

A fairly good description of low-frequency processes of the ion dynamics in the 
shock transition layer can be achieved with hybrid codes (e.g. lWinske fc Omidil ( 19961 ) 
and references therein). Hybrid code modelling, which interprets protons as particles 
and electrons as an inertialess liquid, has made it possible to describe some important 
features of the (sub)sh ock waves at the microscopic scale of some hundred t imes the 
ion inertial length (e.g. lQuestlll988l : iLembeee eraDl2004l : iBurgess et al.ll2005l ). 

A typical initialisation of a shock in the hybrid code simulations is to inject a rela- 
tively cold ion beam (say at the right-hand boundary) and to put a particle reflecting 
wall at the left-hand boundary of a simulation box. In that case the shock is moving, 
and the available simulation time is limited, given the finite size of the box. The limited 
simulation time and the particle statistics per cell are challenging the direct modelling 
of the origin and evolution of the energetic non-thermal particle population in a shock. 
To increase the statistics the m acro-particle splitting method is used (see e.g. lOuestl 
ll992l : lGiacalone fc EUisonlboOOl ). 

In Fig. [2] we show the structure of the magnetic field in a quasi-perpendicular shock 
(inclination angle « 80°) simulated with a hybrid code for the upstream plasma 
parameter /3 ~ 1. The parameter /3 = A4a/A4s, characterises the ratio of the thermal 
and magnetic pressures. The shock is propagating along the x-axis from the left to 
the right. The magnetic field is in the x — z plane. The system is periodic in the y- 
dimension. Phase densities of protons Vx — x, Vy — x, Vz — x are shown in Fig. [3] in the 
reference frame where the particle reflecting wall (at far left) is at rest while the shock 
front is moving. The incoming plasma beam in the simulation was composed of protons 
(90 %), alpha particles (9.9 %) and a dynamically insignificant fraction of oxygen ions 
(Ovil). 

In most of the cases non-relativistic shocks simulated with different hybrid codes 
had the upstream plasma parameter /3 ~ 1. In some cosmol ogical shocks, for e xample in 
hot X-ray clusters, plasma parameter P could be ~ 100 fsee lBvkov et al.ll2008l - Chapter 
8, this volume). The nature of collisionless shocks in the hot low magnetised plasmas 
could be difi'erent from that in case of /3 ~ 1. While the processes of shock formation in 
the high beta plasmas still require a careful study, there are two experimental studies 
establishing the existence of the collisionless shocks for that case. A collisionless shock 
in a laboratory experiment with a laser-produced ablating plasma of /3 ~ 300 was found 
to have a thickness about 10 c/upc, that is orders of magnitude le ss than the Coul omb 
mean free path of both ions and electrons in that system (see iBell et al.|[l988h . In 
space plasma the terrestrial bow shoc k under high beta c onditions was observed with 
the ISEEl and ISEE2 spacecraft by iFarris et all (|l992l '). These measurements were 
compared with and found to be in agreement with the predicted values of the Rankine- 
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Fig. 3 Hybrid simulated proton phase density in a quasi-perpendicular shock (80° inclination). 
The shock is moving from left to right in the reference frame where the particle reflecting wall 
is at rest. The figures show the proton phase densities in Vx — x, Vy — x and Vz — x projections 
from top to bottom respectively. 



Hugoniot relations using the simple adiabatic approximation and a ratio of specific 
heats, gamma, of 5/3. Large magnetic field and density fiuctuations were observed, but 
average downstream plasma conditions well away from the shock were relatively steady, 
near the predicted Rankine-Hugoniot values. The magnetic disturbances persisted well 
downstream and a hot, dense ion beam was detected leaking from the downstream 
region of the shock. 
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3.2 Heating of ions in coUisionless shocks 

The heating processes in coUisionless shocks are non-triviaL The irreversible transfor- 
mation of a part of the kinetic energy of the ordered bulk motion of the upstream flow 
into the energy of the random motions of plasma particles in the downstream flow in 
coUisional shocks is due to Coulomb or atomic particle collisions. In coUisional non- 
radiative shocks without slowly relaxing molecular components, the standard single- 
fluid RH relations are applicable just after a few coUisional lengths. The standard 
single-component shock model predicts a particle temperature kT = {3/l6)mv'^Yi for 
7 = 5/3. 

The particle distributions in the coUisionless shocks are not Maxwellian. Thus, 
instead of the standard equilibrium temperature the appropriate moments of the par- 
ticle distribution function characterising the width of randomised velocity distributions 
are used. Moreover, the particle velocity distributions are typically anisotropic. It is 
clearly seen in hybrid simulated proton phase density: in Fig. [3] the velocity distribu- 
tion widths are different for different projections. One can see in Fig. |3] that a fraction 
of the incoming ions is reflected by the shock magnetic fleld jump providing multiple 
inter-penetrating flows of gyrating ions. Then the field fluctuations randomise the ion 
phases producing a "coarse-grained" distribution characterised by an effective temper- 
ature estimated as the second moment of the velocity distribution. 

An analysis of interplanetary coUisionless shock observations made with Advanced 
Composition Explorer by Korreck et al. ( 2007,) indicated that quasi-perpendicular shocks 
are heating ions more efficiently than quasi-parallel shocks. It was also found that ef- 
fective temperatures of different ions are not necessarily proportional to ion mass, but 
also depend on the shock inclination angle and plasma parameter /3. The widths of the 
coUisionless shocks are extremely narrow (below the astronomical unit) and thus, the 
observed temperatures would depend on the temperature eq uilibration process es (both 
Coulomb and collective) that we will discuss elsewhere (e.g. iBvkov et al.ll20o"^ - Chap- 
ter 8, this volume). The temperature equilibration of different plasma constituents in 
the WHIM can be studied with spatially resolved spectroscopic observations and thus 
is a good test for shock models. 



3.3 Heating of electrons in coUisionless shocks 

Electron kinetics in coUisionless shocks are different from those of ions. Since most of 
the observable emission comes from the electrons, they require a careful study. Shocks 
transfer a fraction of the bulk kinetic energy of the ion flow into large amplitude nonlin- 
ear magnetic fluctuations on a short scale of the transition region (see Fig. [2| . It is im- 
portant that the thermal electron velocities in the ambient medium are higher than the 
shock speed for a shock Mach number A4s < \/mp/mc, allowing for a nearly-isotropic 
angular distribution of the electrons. Non-resonant interactions of these electrons with 
the large- amplitude turbulent fluctuations in the shock tra nsition region could result in 



coUisionless heating and pre-acceleration of the electrons ( Bvkov fc Uvarov|[T999. ). In 
Fig.|4]a simulated electron distribution {p^N(z,p)) is shown as a function of the dimen- 
sionless electron momentum p/ ^2mcTi, where Ti is the initial electron temperature in 
the far upstream {z —oo). The solid curves are the simulated electron distribution 
functions at the left boundary (5 — 0) of the transition region clearly seen in Fig. [2j 
and at the end of the region (z = 1). Dotted lines are the Maxwellian distribution fits 
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Fig. 4 Electron distribution function simulated in the model of electron heating by stro ng ion 
gyroradii scale magnetic fluctuations in a coUisionless shock bv lBvkov fc Uvarovl II1999IV The 
left panel is for a =< {SB/B)"^ >= 0.4, the right panel for a = 0.1. 



allowing to estimate the effective electron temperatures T'^^ measured relative to Ti. 
Note that T^^ = 1.2 at 5 = because of the electron diffusivity effect. One may also 
clearly see the appearance of non-thermal tails indicating a Fermi type acceleration. 
It is worth noting that the presence of large-amplitude waves in the shock transition 
region erodes many of the differences between quasi-parallel and perpendicular shocks, 
making the electron injection mechanism in that model to be similar for these shocks. 

The analysis of observa tional data on both interplanetary and supernova shocks 
bv lChavamian et al.l ( 20071 ') indicates that theelectron heating efficiency i.e. Tc/Ti is a 
declining function of the shock velocity. These authors discussed a model of electron 
heating with a constant level of electr on heating over a wide range of shock speeds (see 
also Fig. 4 in iBvkov fc Uvarovlll999l) . while the ion heating is an increasing function 
of shock velocity . 



3.4 Gas heating and entropy production in weak internal shocks 

Heating/acceleration efficiency with weak and moderate strength MHD shocks can be 
estimated by calculating the energy dissipation rate £h of a- directed gas motion per 
unit area of a weak shock. Defining = Vg\^pTAs, where As is the difference of the 
entropies (per unit mass) behind and ahead of the shock front, one can evaluate the 
thermal gas heating as 



£h = (5/4)(Ms-l)^^sheT, 



(12) 
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where bt is the gas internal energy per unit volume (cf . Eq. Ilip . The energy transfer 
to the reflected nonthermal particles can be estimated from: 

ecR = (-'^s - l)^«sh£B, (13) 



where eb is the magnetic energy density ( Bvkov fc Toptvginlll987l ). From Eg. 1121 and 



Eg. 1131 one can see tha t the gas heating is of third order in {Ms — 1) ^ 1 (cf. 
iLandau fc Lifshitj[l959l ). while the wave damping due to the particle acceleration is 
of second order. Note that in the outer parts of galaxy clusters one would typically 
expect et ^ £b- However, the central regions of such a cluster could have et compa- 
rable to E-Q, as it is the case in the Milky Way. Thus, the weak shocks in the central 
regions could efficiently accelerate nonthermal particles, reducing the heating of the 
gas. Particle acceleration by an ensemble of large scale shocks in a cluster of galaxies 
can create a population of non-thermal particles of sizeable pr essure. This m ay imply 
a non-steady evolution of non-thermal pressure as modelled bv lBvkovl(|200lh . 



4 Energetic particle acceleration in collisionless shocks 

The reflected ions with a gyro-radius exceeding the width of the shock transition region 
can then be efficiently accelerated, via the Fermi mechanism, by converging plasma 
flows carrying magnetic inhomogeneities and MHD waves. In perpendicular shock a net 
transverse particle momentum gain is due to the work of the electric field on the particle 
drift motion. The electric field perpendicular to the shock normal exists in all the 
reference frames for the perpendicular shock. The particle of a momentum p crossing 
back and forth the shock front and being scattered by MHD waves carried with a flow of 
velocity u would undergo a momentum increment Ap « p-^+0{{u/v)'^) per scattering. 
A velocity proflle in the plane shock is illustrated by the dashed line in Fig. [S] in the 
test particle case where one neglects the back reaction effect of accelerated particles 
on the shock. One way to calculate the accelerated particle spectra in a scattering 
medium is to use the kinetic equation in the diffusion approximation. The scattering 
medium in that approach is characterised by momentum dependent particle diffusion 
coefficients ki{p) and k2{p)- The shock is considered as a bulk velocity jump (see the 
dashed line in Fig. [Sjl assuming that the test particles are injected a,t p — pq and the 
gyroradii of the particles are larger than the shock width. Therefore, in a test particle 
case particles must be injected at some super-thermal energy to be accelerated by the 
shock. A solution to the kinetic equation for a nearly isotropic test particle distribution 
in the phase space is a power-law momentum distribution f{p,x) oc {p/po)~^, P ^ Po 
where the index 

b = (14) 
r- 1 ^ ' 



depends on the shock compr ession ratio r jAxford et al.lll977l : lKrvmskiill977l : [Bei]|l978l : 
iBlandford fc Ostrike"il 119781 ). The CR spatial distribution in the model is illustrated in 
Fig. [51 For a strong shock of A4s ^ 1 and A^a 2> 1 the compression ratio given by 
Eqs. [6]and[8]is close to 4 if 7g = 5/3 (or even larger if relativistic gas dominates the 
equation of state). The pressure of the accelerated particles is 

oo 

PcR = -^ / P-" f{p,^)p'^ fip- (15) 
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Then for b — 4 one may see that PcR oc ln(pmax/po) indicating a potentially large 
cosmic ray (CR) pressure, if CRs are accelerated to pmax S> Po- The maximal energy 
of accelerated test particles depends on the diffusion coefficients, bulk velocity and 
scale-size of the system. The finite scale-size of the shock is usually accounted for 
by an energy dependent free escape boundary located either in the upstream or in 
the downstream. For electrons pmax can also be limited by synchrotron (or inverse- 
Compton) losses of relativistic particles. 

The test particle shock acceleration time Ta(p) can be estimated from the equation 

r4p) = ^ /^/fclM^WAdp (^g^ 

Win - U2n J \ "In W2n / P 

Po 

where the normal components of the shock upstream and downstream bulk velocities 
win,'U2„ are measured in the shock rest frame. Estimations based on a more rigorous 
approach distinguishing between the mean acceleration time and the variance does not 
change the results substantially, given the uncertainties in the diffusion model. 



5 Cosmic-ray modified multifluid shocks 

The efficiency of the upstream plasma flow energy conversion into nonthermal particles 
could be high enough providing a hard spectrum of nonthermal particles up to some 
maximal energy e^,. If the efficiency of ram energy transfer to the energetic particles 
is high enough, an extended shock precursor appears due to the incoming plasma flow 
deceleration by the fast particle pressure. The precursor scale L is of the order of 
(c/vsh)-^* ~ orders of magnitude larger than the width of the shock transition region 
(see Fig. [5)). Here A* is the maximal mean free path of a particle in the energy-containing 
part of the spectrum and v^i^ is the shock velocity. We shall later refer to these energetic 
particles as cosmic rays. 

It has been shown that the front of a strong collisionless shock wave consists of 
an extended precursor and a viscous velocity discontinuity (subshock) of a local Mach 
number that is smaller than the total Mach number of the shock wave (see Fig. [5]). The 
compression of matter at the subshock can be much lower than the total compression 
of the medium in the shock wave with allowance for high compression in the precursor. 
We shall refer later to such shocks as CR-modifled. 

The large scale ("macroscopic") structure of a CR- modified shock can be modelled 
by a two-fluid approach with a kinetic description of nonthermal particles (see e.g. 



by a two-tiuid approach witn a kmetic description ot nontnermal particles (see e.g. 
iBlandford fc Eichler|[l987l : iBerezhko et"al]ll99d : iMalkov fc DrurvlboO ll: 'Blasill 2004l and 



references therein) or by a Monte Carlo method (e.g. Joner^r^Uisonl lQQl: Elli soiTet al.l 
In both methods some suitable parameterisation of particle scattering process 
must be postulated a priori. Monte-Carlo simulations, however, have no assumption of 
isotropy for particle distributions, and that allows an internally self-consistent treat- 
ment of thermal particle injection. While the injection depends on the assumptions 
made for the particle pitch-angle scattering, these assumptions are applied equally to 
particles of all energies. The Monte Carlo technique eliminates a free injection param- 
eter, which is present in all models based on the diffusion approximation and is used 
to set the injection efficiency. The strong feedback between injection, shock structure, 
and magnetic field amplification makes this property of the Monte Carlo technique par- 
ticularly important. The Monte Carlo technique allows to iteratively obtain a shock 
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Fig. 5 A sketch illustrating the structure of a cosmic-ray modified shock. The dashed line is 
the shock velocity jump corresponding to the test particle case. The dotted line is a spatial 
distribution of accelerated particles at some momentum p ^ po- The solid line is the CR 
modified shock velocity profile with the precursor and subshock indicated. 



velocity profile and particle distribution function conserving mass, momentum and en- 
ergy fluxes taking into account the nonlinear feedback from the accelerated energetic 
particles. 

In Fig. [6] Monte Carlo simulated proton spectra (multipli ed by [p/{mpc)]'^) are 
shown, in the downstream shock from Iviadimirov et aL I (|2006l 'l. To illustrate the de- 
pendence of the maximal energy of an accelerated proton on the system scale size, 
a free escape boundary condition was applied at some distance from the subshock 
position in the shock rest frame. The heavy solid and dotted curves in the right 
panel correspond to the free escape boundary located at a distance lO^rgi (where 
Tgi — mpUgijC/eBi), the dashed curve has lO^rgi, and the hght solid curve has lO^rgi. 
The simulations were done for a supernova shock in the interstellar medium with a 
shock speed v^h = 5000 kms~^ and an unshocked proton number density ni = 1 cm^'^. 
In the left panel the spectra are given for the same position of the free escape boundary, 
but for different prescriptions of the scattering model. 



5.1 Magnetic field amplification in CR-dominated shocks 

An important predicted feature of strong shocks with efficient CR acceleration is 
the possibility to ampli:^ an initial seed magnetic field by orders of magnitude (e.g. 
iBeU fc LuceMboOll : iBei] |2004| ) . CR current and CR pressure gradient upstream of the 
strong shock could drive magnetic fluctuations on the shock precursor scale length. 
The CR-shock precursor scale L is ~ (c/vsh)-^* which is expected to be above a kpc, 
moreover, the width is L > 100 kpc for a shock of a size comparable to that of a galaxy 
cluster. The precursor scale size L is 3> 10^ times larger than the subshock transition 
region where strong small scale magnetic field fluctuations are directly produced by 
instabilities of super- Alfvenic bulk plasma flows illustrated in Fig. [5] That small scale 
fluctuations are responsible for bulk plasma motion dissipation process and adiabatic 
amplification of the transverse magnetic field in collisionless shocks. At the same time 
the collisionless dissipation process is thought to inject a minor fraction of incoming 
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Fig. 6 Spectra of protons accelerated by a strong shock. The spectra were simulated with a 
non-linear Monte-Carlo model whic h accounts for particle i njection and magnetic field ampli- 
fication by the shock (for details see I Vladimirov et al.|[200^ . On the left panel the transitions 
from the thermal-like component to the high energy tail are marked by pinj . On the right panel 
different curves correspond to different locations of the free escape boundary (see in the text). 



particles to be accelerated to high energies by Fermi mechanism. Recent models of 
diffusive shock acceleration allows a substantial fraction (say, 30 %) of the MIfD shock 
ram pressure to be converted to accelerated particles filling a vicinity of the shock of 
the scale L. The large scale current and density gradient of the accelerated CRs may 
convert a fraction of the CR energy to magnetic field due to multifluid instabilities of 
different kinds providing a way to amplify the initial magnetic field by a factor larger 
than the shock compression ratio. 

Recent non-linear simulations of magnetic field amplificat ion in diffusive shock 
acceleration by a Mo nte- Carlo model I Vladimirov et al] bood ) and a kinetic model 
( Amato fc Blasi 2006) confirmed the possibility of a significant effect. The amplitude of 
the fluctuating magnetic field energy density Wb is of the order of the shock accelerated 
CR pressure which is in turn a substantial fraction of the shock ram pressure 0.5 pi Hgj^. 
Here pi is the shock upstream ambient gas density. 

For typical cluster parameters the discussed mechanism could provide a /iG range 
magnetic field amplitude in a hundred kpc range scale of CR-modified shock precursor. 
The Faraday rotation measure RM provided by a strong CR-dominated shock in a clus- 
ter can reach values of > 10 rad m~^ and even a few times higher. For the case of the 
so-called Bohm diffusion model the rotation measure RM is proportional to the maxi- 
mal energy of the ions in the energy-containing part of the CR-spectrum accelerated by 
the shock. Radio observations, Faraday rotation and synchrotron-C ompton emission 
measurements are used to e stimate the magnetic fields in clusters fe.g. lCarilli fc Tavloij 
I2OO2I : iNewman et"aLl \2Q0i ). Large filaments of polarise d radio emission of scale size 
about 400 kpc were disco vered by Govoni et al.l 1 20051 ) in the halo of the cluster of 



galaxies Abefl 2255 and bv lBagchi et all (|2006l ) in Abell 3376 (see Fig. El). They could 
be connected to large scale shocks due to accretion/merging activity of the cluster. 
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5.2 Gas heating and entropy production in strong CR-modified shocks 

An exact modelling of a collisionless shock structure taking into account the nonthermal 
particle acceleration effect requires the nonperturbative self-consistent description of a 
multi-component and multi-scale system including strong MHD-turbulence dynamics. 
Such a modelling is not feasible at the moment. Instead, a simplified description of a 
multi-fluid strong shock structure can be used with an appropriate parameterisation of 
the extended pre-shock and of the gas subshock. The predicted observable characteris- 
tics of the shocks can be confronted to the observational data. We will now discuss the 
effects of plasma heating by modified shocks and then make some specific predictions 
for possible observational tests. 

In the shocks with efficient high energy particle acceleration the energy flux carried 
away by escaping energetic particles Qosc must be accounted for in the energy continuity 
equations. The energy loss results in a lower effective adiabatic index, but it allows to 
increase the total compression of the gas in the shock downstream. 

The total compression ratio rtot of a strong MHD shock modified by an efficient 
nonthermal particle acceleration can be estimated as 

not = (17) 



7 - \/l + 2(72 - l)Qcsc/pa^;i 



assuming that the energy density in the shock upstream is dominated by the ram pres- 
sure and that the CR escape is through the cut-off momentum regime (e.g. lMalkov fc Drurvl 
I2OO1I ). Here 7 is the effective adiabatic exponent. In Fig.[7]we illustrate the dependence 
of the compression ratio on Qesc/ paVgh ^o'" 7 = 4/3 and 5/3 assuming that the effec- 
tive adiabatic exponent is between the two values depending on the spectrum of the 
accelerated relativistic particles. 

The distribution function of nonthermal particles and the bulk fiow profile in the 
shock upstream region are sensitive to the total compression ratio rtot. Thus, the exact 
calculation of the escape flux Qosc can be performed only in fully nonlinear kinetic 
simulations. Nevertheless, an approximate iterative approach (e.g. in the Monte Carlo 
model discussed above) can be used to make the steady-state distribution function 
consistent with the shock compression assuming some diffusion model. The subshock 
is the standard gas viscous shock of a Mach number A^sub- For that simplified two-Buid 

model of a strong CR-modified shock the effective ion temperature in the downstream 
(2) 

T. can be estimated for a shock of a given velocity, if rtot and rgub are known: 



r « 0(-Ms.b) where ^(A^,.,) = ^/^^^^^'^ 



Single fiuid strong shock heating represents the limit A^sub = Als 3> 1, since there 
is no precursor in that case, resulting in Eq. |9] In single-fluid systems the compression 
ratio rtot = rgub ^ (7g + 1)/ (7g ~ 1) does not depend on the shock velocity and Eq. [18] 
reduces to Eq. [9] However, in multi-fluid shocks the total compression ratio depends 
on the shock velocity and could be substantially higher than that in the single-fluid 
case. This implies somewhat lower postshock ion temperatures for the strong multi- 
fluid shock of the same velocity and could be tested observationally. It is convenient 
to introduce the scaling rtot('Ush) oc v^^ to describe the different cases of strong shock 



heating fsee lBvkovl (|2005l ') for details). Then from Eq.[T8l Ti2 cx (p{Msub) "v^^ The 
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Fig. 7 Total compression ratio r'tot of a strong MHD shock modified by efficient particle accel- 
eration as a function of the energy escape ffux Qesc/^kin carried by energetic particles, where 
^kin = Pl^ah/^' '^'^^ Upper curve (dotted) corresponds to an effective adiabatic exponent 7 = 
4/3 (relativistic gas), while the lower (solid) curve corresponds to 7 = 5/3. 



subshock Mach number Mgub depends, in general, on A^s and Ma.- Thus, an index 
a approximates the velocity dependence of "^(A^sub) oc v^^. Finally, if Ti2 oc w^j^, then 
the index a = 2(l — ^)+a . For t he case of sh ock precursor heating by CR generated 
Alfven waves, the index a 1.25 (|Bvkovll2005h . 



A distinctive feature of multi-fluid shocks is their high gas compression rtotC^sh) 
that could be well above the single fluid shock limit (7g + l)/(7g — 1) (see Fig. [7]). At the 
same time entropy production for a strong multi- fluid shock scales as ''tot('ysh)~^'^'^~'^^^ 
and it is significantly reduced compared to the single-fluid shock of the same velocity. 
The effects are due to energetic particle acceleration and magnetic field amplification. 

Energetic particles penetrate into the shock upstream region. They are coupled 
with the upstream gas through fiuctuating magnetic fields (including the Alfven waves 
generated by the energetic particles). Magnetic field dissipation provides gas preheating 
and entropy production in the extended shock precursor. Such a heated pre-shock 
region of kT < 0.5 keV would appear as an extended filament of width L ~ (c/ugh)-^* ^ 
3 X 10^^e*i3~g cm. Here (in GeV) is the highest energy of the hard branch of the 
accelerated particle spectrum. If B^g ~ 0.1 in the clu ster outskirts and if the hard 
spectrum of energetic nuclei extends to ~ 10^ GeV (cf. iNorman et al]|l995l ) we have 
L ~ 1 Mpc and even wider. Projected on a hot X-ray cluster, such filaments could 
produce a soft X-ray component "excessive" to that produced by the hot cluster. A 
warm gas (~ 0.2 keV) emiss ion filamen t foun d with XMM-Newton in the outskirts of 
the Coma cluster by iFinoguenov et al.l ( 20031 ) could be an extended heated precursor 
of a strong multi-fluid accretio n shock. For a detailed review of the soft X-ray/EUV 
excesses see lOurret et ahlboOSl - Chapter 4, this volume. 
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Fig. 8 Abcll 3376 (taken from Bagchi et al. ( 2006): X-ray emission from XMM-Newton archive 
data, with VLA 1.4 GHz radio contours superimposed. The ellipse shows an elliptical fit to 
the peripheral radio structures, and the '+' marks its centre. The circles mark the positions 
of the two brightest cluster galaxies. 



5.3 ICM entropy production by multifluid accretion shocks 

Cold gas falling into the dark matter (DM) dominated gravitational well passes through 
a strong accretion shock. T he shock is a source of gas entropy production in the inter- 
clust er medium (ICM) (e.g. lKnight fc Ponmanll 19971 : iTozzi fc Normanll200ll : IVoit et all 
I2OO3I ). The post-shock entropy K = K\jT/p^/^ used in the ICM analysis and simu- 
lations (e.g. iBialek et allboOll ) is related to the standard thermodynamic entropy s 
through K oc exp(s/cv). In the standard scenario with a single-fluid accretion shock 
the post-shock entropy scales K^^ oc v^^iPi 

(e.g. lVoit et al]|2003l ). 

The multi-fluid nature of the coUisionless accretion shock modifies the standard 
scaling relation to be 

^mf « «sh[not(^'sh)]"^'+^'=V(^^sub) pI'"''"^- (19) 

The compression ratio in CR-shocks is higher than in a strong single-fluid shock of 
the same velocity resulting in reduced post-shock entropy production. For example, 
in the case of Alfven heating the post-shock entropy of a multi-fluid shock reduces 
as Kj^f/Kgf ~ (15/A1a) for Ala > 15 and A4s > Ma- Here and below in numerical 
estimations we assume jg — 5/3, though a non-thermal baryonic component could 
reduce the index 7g. 
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Since rtot(i'sh) a-nd 0(A1sub) s^re shock velocity dependent, the simple scaling K oc 
^'ihPl ^® '^'^^ valid. In CR-modified shocks K^f oc I's^p^^ or A'mf oc T'^/°- ^ where 
V = 2 — {1 + 7g)^ + cr. For the case of Alfve n wave heating the in dex i/ is < 1.25 and 
Krrtf i s oc T"'^ assuming 7g — 5/3. Recentlv IPonman et al.l ( 20031 ) and IPiffaretti et al.l 
( 20051 ) found that the dispersion in the observed cluster entropy profiles is smaller 
if an empirical r elation K oc is used instead of the standard K oc T (see also 

IPratt et al.ll2006l ;i. 

Consider the simple model of smooth accretion of cold gas through a strong accre- 
tion shock bv lVoit et all (|2003h . The gas of velocity t;ac accretes at a rate Mg through 



the shock at a radius rac where 

Mg = 47rrLpi«ac, «ac = 2GM^r-\ ? = 1 - rac/rta- (20) 
Here M{t) is the cluster mass and rta is the matter turnaround radius. Then the 

(2) 

entropy just behind the multi- fluid shock is expressed through (uac) and 

P2 = rtot{vac)pi- In the Alfven wave heating case 

K^{{t) « {Mt)'-^+'^'^^^, instead of 
Ksf{t) oc (Mt)^/'^ in the single-fluid regime. A multi-fluid shock results in a slower 
post-shock entropy production. As we have noted above, the regime of CR-shock com- 
pression depends on the plasma parameter (3 in the infalling gas. The plasma param- 
eter l3 is currently poorly known because the intercluster magnetic fields are not well 
constrained. The effects of shock modifications are important for both the models of 
smooth accretion of cold gas and for accretion of hierarchical structures. 

Preheating of accreting gas by different physical pr ocesses (e.g. due to ear ly star 
formation in a protocluster region) was suggested by lEvrard fc Henrvl ( 199lh , as a 



possible reason for the breaking of the scaling relations for pure gravitational cluster 
compression bv iKaiseij (|l986l ). The observed high metallicity of clusters at different 



redshifts indicates that strong starburst activity was highly likely at some stage. The 
preheating produce s some initial level of gas entropy ("ent r opy fi oor" , see e.g. exten- 
sive simulations bv iBialek et al.ll200ll : iBorgani etall l200ll . 120051 '). Multi-fiuid strong 
shocks provide a natural alternative way of preheating accreting gas. The non-thermal 
components are essential for detailed modelling of global properties of X-ray clusters, 
including the mass-temperature and luminosity-temperature relations ( Ostriker et al.l 

ioosi). 



6 Shocks in large-scale structure 

Simulations of the cosmic large-scale structure (LSS) predict that about 40—50 % of 
baryons at epoch z < 2 could reside in the Warm-Hot Intergalactic M edium (WHIM) 
with temperatures 10^ — 10^ K at moderate overdensities 5 ~ 10 ( Cen fc Ostriker! 



ll999l : lDave et allboOll ). The WHIM heating is due to shocks d riven by gravitati onally 
accelerated flows in the LSS structure formation scenario (e.g. lKang etaklboOTl ). The 
simulations demonstrate that the X-ray and ultraviolet O vi, O vil and O vill lines and 
the Hi Lya line traces the low-density cosmic web filamentary structures. Interven- 
ing metal absorption systems of highly ionised C, N, O, Ne in the soft X-ray spectra 
of bright AGN were suggested to trace the WHIM. The detection of shocked WHIM 
requires very sensi tive UV and X-ray detectors, both for ab sorption and for emission 
processes (see e.g. iKaastra et al.l [2OO8I : iRichter et al.l I2OO8I - Chapters 9 and 3, this 
volume). Dedicated future missions like the Cosmic Origin Spectrograph (COS), the 
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Fig. 9 Metallicity map of the intracluster gas in Abcll 3376 (taken from lBagchi et al 



X-Ray Evolving Universe Spectrometer {XEUS), Constellation-X and the Diffuse In- 
tergalactic Oxygen Surveyor (DIOS) will provide high resolution spectroscopy of the 
shocked WHIM. The WHIM ions of different charge states have highly non-equilibrium 
(anisotropic) initial states just behind a collisionless shock that relaxes to equilibrium 
states through Coulomb collisions. As was discussed above a strong collisionless shock 
could generate a spectrum of MHD-fluctuations. These MHD-fluctuations can carry 
a substantial fraction of the shock ram pressure. The velocity fluctuations will result 
in non-thermal broadening of the lines, potentially important for simulations of emis- 
sion/absorption spectra of the WHIM and observational data analy sis. Specific feature s 
of collisionless shock heating of the WHIM ions are discussed in iBvkov et al.l |2008| - 
Chapter 8, this volume. In this paper we discuss only a few observations of clusters of 
galaxies. 



6.1 Evidence for shocks in galaxy clusters 

Clusters of galaxies are believed to form within the hierarchical build up of the large 
scale structure of the Universe. Small objects collapse first and then merge in a complex 
manner to form larger and larger structures. Therefore, once in a while during their 
formation, clusters of galaxies undergo so called major merger events. In such events, 
proto-cluster structures of similar masses (typical ratios 1:10 — 1:1) are colliding with 
super sonic velocities (typically several 1000 kms~^). These merging events are a source 
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of shocks and turbulence. They redistribute and amplify magnetic fields, and they 
are a source of acceleration of relativistic particles within the intracluster medium. 
Additionally, it is expected that the accretion of the diffuse, unprocessed (and therefore 
relatively cold) matter onto the DM node of the cosmic web creates a virialisation shock 
(also called accretion shock), which is expected to be located far in the cluster periphery 
(typically a few Mpc from the cluster centre for massive systems). Some examples of 
relevant cosmological simulations are presented bv lOolag et al.|[2008l - Chapter 15, this 
volume. 

X-ray observations, revealing the thermodynamical state of the intracluster medium 
are therefore the natural means for searching for the signatures of such non thermal 
phenomena. However, due to biases in the observational processes, caused by the com- 
plex temperature structure of the intracluster medium, such si g natur es are very hard 
to detect (for a more detailed discussion see iMazzotta et al. I (|2004l ') and references 
therein). Nevertheless, some detections of shocks in galaxy clusters have been revealed 
by high resolution Chandra an d XMM-Newton observations. For a recent review see 
iMarkevitch fc VikhhninI l|200'if ). One of t he most spectacular exa mples of a merging 
galaxy clusters is the case of IE 0657-56 (|Markevitch et al.ll2002l ). 

The sensitivities of current X-ray instruments are not sufficient to map the state 
and structure of the intracluster medium in the periphery of galaxy clusters. However, 
the discovery of arc like radio emission in the periphery of some clusters (so called 
radio relics), are thought to trace shocks running through the intracluster medium. 
Spectacular examples are Abell 3667 (see also ,Fcrrari ct al. 200 8 - Chapter 6, th is 
volume) or the radio relics in Abell 3376 recently discovered bv lBagchi et al.l ( 20061 ) ■ 
see Fig.[8l Thereby, such radio observations are currently the only possibility to observe 
shocks outside the central regions of galaxy clusters . A more detailed discussion of 
numerical models of such radio relics is provided by bolag et al.ll2008l - Chapter 15, 
this volume. 

Shock s are also exp ected to trigger star formation, as indicated by numerical simula- 
tions (e.g.|Bekkii ri999l ). which will leave detectable imprints in the intracluster medium 
even long (several Gyr) after the shock passed through. Multiple supernova explosions 
in the star forming regions (superbubbles) will ad ditionally pro duce copious small scale 
shocks and accelerate non-thermal particles (e.g. lBvkovll200lh . One of the tracers for 
these processes can be excess metalli city in the intracluste r medium, produced by the 
enhanced star formation period (e.g. ISchindler et al.ll2005l ). Fig. |9] shows the inferred 
metallicity map fo r Abell 3376, indicating previous merger activity of the cluster (see 
iBagchi et al]|2006l ). 



7 Summary 

Cosmological shocks convert a fraction of the energy of gravitationally accelerated ffows 
to internal energy of the gas. They heat and compress the gas and can also accelerate 
energetic non-thermal particles and amplify magnetic fields. We discussed some specific 
features of cosmological shocks. 

• The standard Rankine-Hugoniot relations based on the conservation laws for a 
steady single-fluid MHD shock allow to calculate the state of the fluid behind the shock 
once the upstream state and the shock strength are known. The coplanarity theorem 
for a plane ideal MHD shock states that the upstream and downstream bulk velocities, 
magnetic flelds and the shock normal all lie in the same plane. 
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• Cosmological plasma shocks are likely to be coUisionless as many other astro- 
physical shocks observed in the heliosphere and in supernova remnants. We review the 
basic plasma processes responsible for the microscopic structure of coUisionless shocks. 

• CoUisionless shock heating of ions results in a non-equilibrium state just behind 
a very thin magnetic ramp region with a strongly anisotropic quasi-Maxwellian ion 
distributions. The possibility of coUisionless heating of electrons by electromagnetic 
fluctuations in the magnetic ramp region depends on the extension of the Uuctuation 
spectra to the electron gyro-scales, and could depend on the shock Mach number. Then 
the Coulomb equilibration processes are operating on the scales much larger than the 
coUisionless shock width. 

• Extended MHD shock waves propagating in turbulent media could accelerate 
energetic particles both by Fermi type acceleration in converging plasma flows and by 
DC electric field in quasi-perpendicular shocks. If the acceleration is efficient, then the 
strong shock could convert a substantial fraction (more than 10 %) of the power dissi- 
pated by the upstream bulk fiow to energetic particles (cosmic rays) . The compression 
ratio rtot at such a shock can be much higher, while the ion temperature behind the 
shock oc and the post-shock entropy are lower, than that in a standard single fiuid 
shock. The shock structure consists of an extended precursor and a viscous velocity 
jump (subshock) indicated in Fig. [5l 

• Strong coUisionless plasma shocks with an efficient Fermi acceleration of energetic 
particles could generate strong MHD waves in the upstream and downstream regions 
and strongly amplify the upstream magnetic fields. A distinctive feature of the shock 
is a predicted possibility of gas pre-heating in the far upstream region due to MHD 
wave dissipation, that can produce an extended filament of temperature > 0.1 keV. 

• Shock waves both from the cosmic web formation processes and those due to 
cluster merging activity can play an important role in clusters of galaxies. Direct evi- 
dences for such shocks, as traced by radio relics and the temperature jumps in X-ray 
observations havebeen found only in a small number of clusters, and thus we need more 
observations. 
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